Simulations of multicomponent relativistic thermalization

Multicomponent relativistic fluids have been studied for decades. However, simulating the dynamics of particles and in such a mixture has challenge due to fact that simulations are computationally expensive three spatial dimensions. Here, we report on development application multi dimensional Monte Carlo code explore thermalization process multicomponent environment inexpensive way. As an illustration simulate fully Brownian-motion-like solution high mass particle (proton) bath low-mass (electrons). We follow ultimate equilibrium distribution Brownian-like as can happen cosmic plasma during big bang nucleosynthesis. also energetic injected into occur, example, by decay massive unstable bang.